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200-0032 

METHOD FOR RECONSTRUCTING 
THE TOPOLOGY OF A POLYGONAL SOUP 

Technical Field 

The present invention relates to 
5 visualization tools for interactive simulations with 
complex vehicle models. 

Background Of The Invention 

To meet the need for interactive simulations 
with complex vehicle models, a generation of 

10 visualization tools have been developed. However, most 
of these tools have limitations inherent in their data 
structure. In this regard, most visualization tools 
have adopted a polygonal soup, that is, an unstructured 
collection of polygons, as their standard data 

15 representation. However, a variety of existing 

applications, such as CFD {Computational Fluid 
Dynamics) , FEA (Finite Element Analysis) , and assembly 
simulation, require more than a cloud of points or 
triangles . 

2 0 Thus, a need exists for a method of 

reconstructing the topographical information 
(relationships between polygons) of a polygonal soup in 
order to overcome the known problems with visualization 
tools . 



2 5 Summary of the Invention 

The object of the present invention is to 
provide an improved visualization tool for interactive 



simulation with complex vehicle models. It is another 
object of the present invention to provide a 
visualization tool that comprises more than a cloud of 
points or tr±angles . 
5 - - In accordance with the present invention, an 

algorithm is provided which automatically reconstructs 
topological information for a given mesh and then 
alters the mesh by introducing, deleting, or splitting 
existing polygons when needed. 

10 In accordance with the present invention, a 

new OctTree space decomposition is used to achieve a 
log2 -complexity search which locates the closest vertex 
in the polygonal soup to a given point in space . A 
technique with linear complexity is then used to locate 

15 all of the triangles connected to that vertex. A 
technique with linear complexity is also used to find 
all triangles connected to a given triangle. The 
triangles are then split to enforce conductivity. 

These and other objects, features, and 

2 0 advantages of the present invention will become 
apparent from the following detailed description of the 
invention when viewed in accordance with the 
accompanying drawings and appended claims. 



Brief Description Of The Drawings 

2 5 FIGURE 1 depicts an exemplary cube 

representation for modeling tessellated meshes; 

FIGURE 2 sets forth index representations of 
the triangles and vertices of the cube depicted in 
Figure 1 ; 



2 



FIGURE 3 depicts a dynamic vector vertex- 
neighbor table; 

FIGURE 4 depicts tables generated from the 
vertex-neighbor table; 
5 - ' FIGURE 5 depicts a dynamic vector edge- 

neighbor table; 

FIGURE 6 illustrates an object compound of 
unconnected strips of connected triangles; 

FIGURE 7 depicts a generalized example of a 
10 step in synchronizing triangle strips; 

FIGURE 8 depicts the formation of new 

triangles ; 

FIGURE 9 illustrates the exemplary cube as 

modified; 

15 FIGURE 10 illustrates another step in the 

sorting process; 

FIGURE 11 illustrates another step in the 
sorting process; 

FIGURE 12 depicts another sorting example; 
2 0 FIGURE 13 depicts still another sorting 

example ; 

FIGURE 14 illustrates an OctTree example; and 
FIGURE 15 depicts a further sorting example. 



Description Of The Preferred Embodimeiit(s) 

25 As indicated, there is a need for interactive 

simulation with complex systems, such as airplanes and 
cars, and in fact visualization tools have been 
developed which are capable of rendering large models 
in real-time. However, in developing the 3-D rendering 
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tools, insufficient attention has been paid to 
architecting systems capable of evolving into effective 
modeling tools. One of the critical limitations is 
inherent in the data structure. Almost all 

5 visualization tools have adopted polygonal soups as 
their standard data representation. However, a variety 
of existing applications, ranging from CFD, FEA and 
assembling simulations, require more than a cloud of 
points or triangles. 

10 Currently, different data representations are 

used to model tessellated meshes. Figure 1 illustrates 
one example represented by a cube 10. The cube is 
modeled as a set of eight triangles, two per side, and 
eight vertices, Vq through . 

15 In the STL (stereo lithography) format, the 

model is described in either ASCII or binary format, as 
a list of triangular elements. The actual coordinates 
of the three vertices of each of the triangles are 
expressly stored and are shown in the lefthand table 

2 0 set forth in Figure 2. The actual coordinates of the 

three vertices are provided on the left side of the 
table in Figure 2 , and the cube is represented as a 
list of 3 X 12 vertices. In this regard, other 
representations, such as Virtual Reality Modeling 
25 Language (VRML) , OpenGL, Optimizer, DirectX, and 
DirectModel, offer other efficient, indexed 
representations . 

Figure 2 shows two separate lists, a vertex 
list with the coordinates of all vertices, and a 

3 0 triangle list where the indexes of the vertices for 

each triangle, relative to the vertex list, are stored. 
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The STL representation is, in fact, a true 
polygonal soup, that is, the geometry is a set of 
unconnected triangles. The index representation is 
better. In this regard, the same vertex is shared 
5 among -multiple triangles and provides some topological 
information. However, the facets are usually connected 
in strips or fans for graphics performance, without a 
collective topological map and with multiple 
representations of the same vertex on different strips. 

10 This causes small gaps in the geometry to appear 
because of numerical approximations. 

Real-time tools, such as collision detection 
tracking algorithms, take advantage of spatial 
coherency between successive time samples. This 

15 requires information about vertex and edge 
connectivity, such as the neighbors of a given triangle 
and of a given vertex. For example, in the case of the 
cube shown in Figure 1, the triangle T4 is connected to 
triangles T^, ^^^d T5, and the vertex V2 is shared 

2 0 among the triangles T-^, T^ , and T^ . 

The present inventive method basically 
consists of three steps or phases. In the first phase, 
the vertex and edge conductivity information is 
developed starting from available information. The 

2 5 second and third phases deal with imperfect meshes. 

During the second phase, vertex duplicates are 
discovered and eliminated. During the third phase, the 
model is re-meshed to realign strips of triangles that 
do not share common vertices. 

3 0 The initial phase is divided into three 

additional steps. In the first step, the index 
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representation of Figure 2 is generated. An initial 
attempt of reducing duplicated vertices is performed at 
this stage. In this regard, the first step operation 
is only required when the geometry is described as an 
5 explicit set of triangles, such as STL. In the second 
and third steps, the vertex-neighbors, and edge- 
neighbor tables are developed. 

The vertex-neighbors table is a dynamic 
vector where, for each vertex, the list of all the 

10 connected triangles is immediately available, as is 
shown in Figure 3 . The algorithm sequentially browses 
the triangle list. For each triangle, the table listed 
in Figure 2 contains the indexes of the three vertices. 
An identifier, pointing to the triangle sequential 

15 number in the triangle-list is added to the index- 
neighbor list, in correspondence to each vertex. For 
example, in the case of the cube shown in Figure 1, 
this method would start by looking at the triangle Tq 
with vertices V^, V,^, and V3 . Table A in Figure 4 is 

2 0 then generated. 

Proceeding with Tj^, the additional entries are 
added to Table B. Thereafter, the complete vertex- 
neighbor Table C is developed. Since the triangle- 
index table is traversed only one time, this creates an 
25 efficient, linear-complexity algorithm. In this 

regard, the number of computations is n^ x 3 , where n^ is 
the number of triangles. 

The edge -neighbors table is a vector of lists 
with an entry for each triangle. Each list has three 

3 0 elements, one for each edge of the given triangle, Vj^;^- 

Vi2, V;^j-Vj^3, and 'V^^-V^^, respectively, depending on the 
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position of the element in the list. These elements 
can be equal to either -1, to represent an unconnected 
edge, or to the sequential number of the connected 
triangle to the selected edge, as can be seen in Figure 
5 5. -Thereafter, the algorithm sequentially browses the 
triangle list, starting with the first element, and 
determines if there are other triangles that share any 
of its three edges. For example, with the cube 
illustrated in Figure 1, the first edge VoV^ of the 
10 first triangle Tq can be viewed as follows: 



Tri angl e 




Vl2 




To 




Vl 


V3 



Two triangles are connected at an edge if 
they share the same two vertices making up the edge 
15 extremes, in this case Vg and . Then, a list of 
potential candidates is found looking for all triangles 
connected to Vg that are not Tg in the vertex-neighbors 
table. This is shown as follows: 



Vertex 












Vo 


To 


Ta 


T9 


Tio 


Til 



20 

Since there are now four triangles, Tg-Tj^^, it 
is necessary to look at the triangle table to see if 
any of the candidates contains the vertex as well: 
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Triangle 


^11 




















Vo 


Tio 


V3 




Vo 


"^11 


V4 




V7 



(3) 



In this particular example, the first edges 
of To and the first edge of Tg share the same vertices 
and Vg . The two triangles are connected, and the 
5 edge -neighbor table is then updated as shown below: 



Triangle 


Til 


^i2 




To 










To 







As a result, after browsing all of the edges 
of all of the triangles, the following table is 
10 completed: 



Triangle 


Til 


Ti2 




To 


Tg 


Ti 


Tio 


Ti 


T2 


To 


T4 




T4 


T3 


Ti 


T3 


Tio 


T2 


T7 


T4 


Ti 


Ts 


T2 


T5 


T7 


T4 


Tg 


Tg 


Til 


T7 


T9 


T7 


Ts 


Tg 


T3 


Te 


To 


Tg 


T5 




Ts 


Tg 


Til 




T3 


Til 


To 


Til 


T9 


Tio 


Ts 


As shown in 


Table (5) 


it is 


noted that 



are not any unconnected edges and, therefore, this 
5 particular geometry is a manifold body. 

The total number of computations depends on 
the average number of triangles connected to a given 
vertex. The triangle/vertex ratio n^/^ ranges from a 
value of three in the case of a tetrahedron, to 4.5 for 
10 the cube described above, to about 5.0 for typical 
automotive models. This means that the number of 
computations to find the triangle connected to a given 
edge is: 

(n,/,-l)-3 (6) 
15 Since each triangle has three edges, and there are n^ 
triangles, then the total number of computations is: 

n^- (nt/,-1) -9 « nt-36 (7) 
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This provides an efficient, linear-complexity 
algorithm. 

The algorithms used to reconstruct the 
topographical information of a tessellated mesh, as 
5 discussed above, rely on two assumptions. The first 
assumption is that duplicate vertices have the same 
coordinates. This requirement is a foundation of a 
variety of search engines, usually based on balanced 
trees or skip-lists, used to detect duplicate vertices. 

10 The second assumption is that the tessellation is built 
without gaps. This translates in having strips or fans 
of triangles all topologically connected at the same 
set of vertices. However, in particular with VRML and 
OpenGL files, the objects are composed of unconnected 

15 strips of connected triangles, as shown in Figure 6. 
This prevents building of edge connectivity since two 
triangles are considered connected at an edge only if 
they share the entire edge. 

As a result, existing visualization tools for 

2 0 complex systems have a plurality of gaps, poor 

alignment, and degenerate meshes. Also, duplicate 
vertices and the synchronization of triangle strips are 
not sufficiently taken into account. 

In the present invention, these problems with 
25 existing systems are considered and taken into account. 
Duplicate vertices are eliminated and triangle strips 
are synchronized. The algorithm for accomplishing this 
is described below. 

If the vertex is duplicated and the copies 

3 0 replaced by the vertex V^^ in the cube model described 
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above to describe the triangles T^^ and T^.^, the new 
vertex and triangle list are then modified as follows: 



Vertex 


X 


Y 


Z 


Triangle 


Vii 


Vi2 Vi3 


Vo- - 


-1 


-1 


+ 1 








Vob 


-1 


-1 


+ 1 


Tio 


V3 


V7 Vob 










Til 


V4 


Vob V, 



(8) 



Due to the new duplicate vertex, the vertex-neighbors 
and edge-neighbors tables are also changed to the 
following : 



Triangle 


Til 


Ti2 


Ti3 


Vertex 


Til 


Ti2 ... . 




To 


Ts 


Ti 


-1 


Vo 


To 


Ta T9 












Vob 


Tio 


Til 




T9 


Ts 


Ta 


-1 










Tio 


T3 


Til 


-1 


J-ii 


-1 


Tio 


Ts 



As a result, all four triangles, To_ Tg, T^q, and T^^, now 
each have an unconnected edge. 

In order to remove duplicate vertices, the 
algorithm first browses the edge table searching for 
unconnected edges. Once an edge is found, duplicate 
copies of either one of the two vertices at the ends of 
the segment are looked for, by searching in the vertex 
list for the closest vertex. In Table 8 above, 
vertices Vg and V^^, have the same coordinates, and 
therefore have minimum distance, equal to zero. 
However, in other cases, the distance is small, but 
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significant. It is necessary to make sure that the mesh 
will not degenerate once a vertex is replaced with 
another one . 

Once a duplicate vertex has been found, the 
procedure is as follows: (a) in the triangle index 
table, all references to the duplicate vertices are 
replaced with a reference to the original vertex; (b) 
the entry in the vertex table for the duplicate vertex 
is removed and the appropriate connected triangles are 
added to the entry list of the original vertex; and; 
and (c) the edge table for all triangles connected to 
"original" vertices are rebuilt. When applied to the 
copy of duplicated vertices, and V^^, this sequence 
restores the original vertex, triangle, vertex- 
neighbors and edge -neighbors tables. 

The problem with synchronizing triangle 
strips, such as those depicted in Figure 6, can be 
generalized as shown in Figure 7 as a vertex falling 
inside another triangle edge. The introduction of two 
new triangles, namely T^^ and T^i, as well as a new vertex 
Vg, changes the vertex and triangle table as follows: 



Triangle 


Vil 


V,2 








Va 




Tib 


V2 


V3 





As evident from the new edge-neighbor table, there are 
three edges that are not connected: 



Tri angl e 






T.3 


To 


Ts 


-1 


Tio 




Ti, 


-1 


T4 


Tib 


T2 


-1 


Tia 



There are different alternatives for building 
a topological connection. First, an artificial vertex 
could be introduced in the triangle Tq . However, this 
would require dealing with polygons having different 
numbers of sides since with a new vertex, T^, would be 
a quadrilateral . Secondly, the vertex Vg could be moved 
to overlap either one of the two vertices or V3, and 
then be eliminated. However, this approach would only 
modify the geometry, with the potential loss of 
important details, if the triangles connected to Vg do 
not share the same normal. Thirdly, the triangle T^ 
could be split into two new triangles, replacing the 
edge VjV^ with V^Vg and VjVg, respectively. 

The present invention implements the latter 
alternative. As in the method dealing with duplicated 
edges, the algorithm browses the edge table looking for 
unconnected edges. Once an edge is found, it then 
looks for any vertex that falls very close to the edge 
itself, but not on one of its vertices. In the example 
shown in Figure 7, the first unconnected edge that is 



Vi2 Vu 



Triangle 




Ti2 


Ti3 


To 


Ta 


-1 


Tio 




Tib 


-1 


T4 


Tib 


T2 


-1 


Tio 
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encountered is the second edge Vj^Vj of the triangle T^: 
Searching in the vertex list, the closest vertex to 
this segment is Vg. Since the distance to the segment 
is equal to zero, and the distance to both segment 
5 extremes" is equal to Vs , the vertex is recognized to 
fall on an inside edge. Once the closest vertex is 
found, it is necessary to search in the vertex- 
neighbors table to see if there is a triangle connected 
to it which has an unconnected edge parallel to the 
10 original edge. Moreover, the triangle must not overlap 
the original triangle. In the "cube" example discussed 
above, there are two triangles connected to Vg : 



Vertex 


T,, T,, 







Once a triangle, edge, and vertex have been 
identified, the geometry is then re-meshed by splitting 
the triangle. First, the triangle table is updated to 
reflect the fact that the original triangle is split 
into two new elements. In general, the original 
triangle is kept with one of its vertex indexes 
replaced and a new triangle is added. The triangle 
table will then change as follows: 





V,2. Vl2 .^^»..yi = 


Vxl Vi2 Vu 


To 


Vo Vi V3 ^0 


Vo Vi Vs 
Vo Vs V3 


il2 



(14) 
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In the generic case, the algorithm must take 
into account which edge the unconnected vertex falls 
on, as shown in Figure 8. The generic rules to update 
the triangle table are then summarized as follows: 







Before 


After 




Triangle 


Vii 


Vi2 




Vii 


Vi2 


Vi3 


edge 

Vu-Vi2 


Tnew 




Vi2 


Vi3 


Vii 

Vnew 


Vnew 
Vi2 


Vi3 
Vi3 


edge 

Vi2-Vi3 


Torig 
Tnew 


Vii 


Vi2 


Vi3 


Vii 
Vii 


Vi2 
Vnew 


Vnew 
Vi3 


edge 

Vi3-Vu 


Tnew 


Vii 


Vi2 


Vi3 


Vnew 
Vii 


Vi2 
Vi2 


Vi3 
Vnew 



(15) 



Thereafter, the vert ex -neighbor table is 
15 updated. In the case of the cube shown in Figure 7, 
this would change as follows: 



25 





Til 


Ti2 








Vo 


To 


Ts 


T9 


Tio 


Til 


Vi 


To 


Tia 


T4 


Ts 


Ta 


Vj 


To 


Tib 


T2 


T3 


Tio 


Va 


Tia 


Tib 










Vertex 


Til 


Ti2 








Vo 


To 


Ta 


T9 


Tio 


Til T12 


Vi 


To 


Tia 


T4 


Ts 


Ts 


V3 


Tib 


T2 


T3 


Tio 


T12 


Vs 


Tia 


Tib 


To 


T12 





(16) 
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In the general case, the following table 
summarizes the rules needed for determining the new 
table, depending on which of the three angles is 
affected by the new vertex: 







After 




Vertex 




Til 


edge 


Vii 








Vil-Vi2 


Vi2 


Torig 








Vi3 




Tnew 






Vnew 






Tnew 


edge 


Vii 








Vi.2-Vi3 


Vi2 










Vi3 


Torig 


Tnew 










... Torig 


Tnew 


edge 


Vii 


Torig - 


Tnew 




Vi3-Vii 


Vi2 




Tnew 






Vi3 










Vnew 




Torig 


Tnew 



(17) 



The last step is to update the edge -neighbors 
table in correspondence of the three vertices of the 
original triangle V^^, ^.s well as the vertex 

projecting in the middle of the edge V^^^, using the 
algorithm described above in equations (6) and (7) . In 
the case of the present example, the edge -neighbors 
table is updated as follows: 



Tri angl e 

Tib 



T4 



(18) 
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The geometry is now manifold as shown in Figure 9. 

An efficient closest -vertex search is 
utilized in accordance with the present invention. If 
the entire geometry is searched, the number of 
5 calculations required to eliminate duplicate vertices 
and split triangles in order to obtain a connected, 
closed manifold is a function of the square of the 
number of vertices and shown by the following equation: 

a _ « «e 2 _ "-e/v 2 _ <ar 3 n^^ 2 

-•"e-«v - 2-~-"v "2"^"'"" ~2"2~2~'"'' ^^^^ 

10 

where ng is the number of edges, r\ is the number of 

vertices, is the edge/ vertex ratio, = — — is the 

2 2 3 2 

triangle/vertex ratio, and a is the percentage of edges 
that are not connected and require some action. Since 

15 a typical tessellated model of a single vehicle 
component usually consists of hundreds of thousands of 
triangles, the number of computations quickly grows and 
becomes unmanageable when dealing with complex 
assemblies. Different solutions are available. Most 

2 0 of these are based on the same philosophy where all of 
the elements of the model are assigned to some bounding 
volumes of simple geometry. In order to reduce the 
number of computations from a linear to a logarithmic 
dimension, the hierarchy of bounding volumes is often 

2 5 used. Some implementations subdivide the space in 
spherical - containers , while others make use of 
variable-side bounding boxes. 

With the present invention, the space is 
decomposed using an OctTree structure. This is a data- 
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structure, similar to a voxel map, and is in current 
use in computer -graphics . An example illustrates its 
concept. If the parent mode of the OctTree is a box 
enclosing the entire workspace, vertices are assigned 
5 to -it.- As soon as there are more than a constant, 
assigned number of vertices in the parent node, for 
example 4, the box is split into a set of 8 identical 
smaller boxes, where the intersection is the empty set 
and the union is the original box. The elements 

10 previously assigned to the original box are then 
reassigned to the smaller containers . This is shown in 
Figure 10. As the sorting continues and new vertices 
are assigned to the proper boxes, the sub-box that is 
filled is once more subdivided into eight smaller 

15 volumes and so on. This is shown in Figure 11. 



node has a pointer linked to either a list of vertices 
or a list of sub-boxes. On the other hand, each vertex 
has embedded in its coordinate the address of the box 

20 containing it. Since the boxes are always split and 
each side is dissected into two equal segments, it is 
sufficient to find the address of the box to normalize 
the coordinates by the size of the workspace and then 
multiple the result by 2", where n is the maximum number 

2 5 of layers, and round off to the nearest integer. For 
example, given the point shown in Figure 12, the 
address in an n=4 layers OctTree is as follows: 



These indexes can be used to find the 
30 containing box. If the parent node points to a list of 
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The data structure is bi-directional. 



Each 




(20) 



vertices, then the search is completed. Otherwise, the 
Most Significant Bits (MSB) defines the address of the 
sub-box which in this case is box 0,0. If this box had 
been subdivided, the child container is found by 
5 shi-fting the indexes to the left and using the MSBs 
once more. In the particular example shown above, the 
procedure ends after the indexes had been shifted three 
times. The sequence of boxes is then as follows: GO, 
00, 10, 11. 

10 Most vertices can be eliminated by looking 

only at those inside boxes whose intersections with the 
given shape is not null. To identify these boxes, one 
possibility is to find the smallest box in the OctTree 
index space enclosing the edge plus a neighborhood 

15 region function of epsilon. To do so, it is first 
necessary to map the edge equation from Cartesian to 
the OctTree index space. The edge is defined as the 
parameter ex as: 

(Xo,yo,Z(j) + a- (Ax,Ay,Az) , with j| Ax, Ay , Az || 2 = 1 
20 and 0 < a < and it starts in Pq = (Xo,yo,Zo) (21) 
and ends in Pj. = (Xo,yo,Zo)+ {Ax,Ay,Az) . 

Assuming all positive Ax,Ay,Az', the maximum and minimum 
indexes of a box, aligned with the system axis and 
enclosing the edges are: 

25 




(22) 
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Since the point must fall in an epsilon- 
neighborhood of the segment, the delta is calculated in 
the OctTree index space as follows: 



(23) 



and then subtracted and added to the edge-enclosing 
box. The closest elements will be searched among all 
the vertices belonging to any of the sub-boxes 
10 satisfying the requisite equation: 

T^^ - {Ai, Ai, Ai,) < T < 4^ + ( Az, Ai, Ai) (24) 
In the example illustrated in Figure 13, this 
would require checking all boxes whose indexes satisfy 
the following equation. 
15 0011,0001-l,l<r<0111,0101 + l,l=>0010,0000<r <1000,0110 (25) 
The sequence of boxes to be checked is 
generated. It is sufficient to traverse the tree, 
depth-first. For the particular example described 
above, the tree is shown in Figure 14. The search 
20 starts browsing the leftmost branch of the 'tree until 
all leaves are traversed, and then continues exploring 
the second leftmost branch of the tree. The last two 
branches are discarded because their indexes are 
outside the allowed range, that is, the following 
25 conditions are not satisfied: 

0010,0000<lxxx,0xxxx<1000,0110 and (26) 

0010,0000<lxxx,lxxxx<1000,0110 (27) 



20 



A more efficient alternative is to only check 
the boxes that are in the close neighborhood of the 
edge, as shown in Figure 15. The equation of the edge 
in Cartesian space, which is 

P(a)= (Xt5,yo,Zo) + a- (Ax,Ay,Az) 

(28) 

with ||Ax, Ay, Az||2 = 1 and O^a^a^ 
becomes the OctTree index space 
T(a) = 



(29) 



10 Using this equation, it is now possible to 

traverse the tree and for each box to determine if it 
contains a portion of the edge. Even if the technique 
requires additional computations to determine the list 
of boxes to be included in the search, it culls some of 

15 the boxes that would otherwise have been selected and 
as a consequence reduces the number of vertices whose 
distance from the target needs to be evaluated. On 
average, the overall number of calculations is less 
than in the previous implementation. Since OctTree are 

2 0 three-dimensional binary trees, the searching algorithm 
in accordance with the present invention has a 
logarithmic, base 2, complexity. 

While the invention has been described in 
connection with one or more embodiments, it is to be 

2 5 understood that the specific mechanisms and techniques 
which have been described are merely illustrative of 
the principles of the invention. Numerous 
modifications may be made .to the methods and apparatus 
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described without departing from the spirit and scope 
of the invention as defined by the appended claims. 
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What Is Claimed Is: 



1. A method for reconstructing topological 
information for a mesh, said mesh comprising a 
polygonal soup of triangles with sides and vertices, 
said method comprising the steps of: 

building vertex and edge connectivity data; 
finding duplicates of vertices; 
removing said duplicates of vertices; and 
realigning strips of triangles without common 

vertices . 

2 . The method as set forth in claim 1 
wherein said step of building vertex and edge 
connectivity data comprises the steps of: 

generating a representative index; 

creating a vertex-neighbor table; and 

building an edge-neighbor table. 

3 . The method as recited in claim 2 wherein 
said step of generating a representative index 
comprises eliminating at least one duplication of 
vertices . 

4 . The method as recited in claim 2 wherein 
said step of removing duplicate vertices comprises: 

searching for unconnected sides of triangles; 

searching for duplicates of the vertices at 
the ends of said unconnected sides; 

replacing all duplicate vertices with 
original vertices; 

adding triangles connected to the duplicate 
vertices to said original vertices; and 
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rebuilding said edge-neighbor table for all 
triangles connected to said original vertices. 

5. The method as set forth in claim 4 
wherein said step of adding triangles comprises 

5 splitting the triangles into new smaller triangles. 

6. The method as set forth in claim 4 
wherein said step of searching for duplicates comprises 
searching in said vertex-neighbor table for the closest 
vertex. 

10 7. The method as set forth in claim 6 

wherein said step of searching for the closest vertex 
comprises using an OctTree structure. 

8. The method as set forth in claim 6 
wherein said step of searching for the closest vertex 

15 comprises using a log2 -complexity search method. 

9 . The method as set forth in claim 8 
wherein said log2 -complexity search method comprises 
using an OctTree structure. 
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Abstract Of The Disclosure 

A method for automatically reconstructing 
topographical information for a given mesh, altering 
the mesh by introducing, deleting, or splitting 
5 existing polygons when needed. An OctTree space 
decomposition is used to achieve a log2 -complexity 
search method to find the closest vertex in the 
polygonal soup to a given point in space. Linear 
complexities are used to find triangles connected to a 
10 given vertex and all triangles connected to a given 
triangle. The triangles are split to enforce 

conductivity . 
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solid CUBE 
facet 

outer loop 
vertex -1 -1 +1 
vertex +1-1+1 
vertex -1 +1 +1 
endloop 

endfacet 

facet 
outerloop 
vertex +1 +1 +1 
vertex -1+1+1 
vertex +1-1+1 
endloop 

endfacet 

facet 

outer loop 
vertex +1 +1 +1 
vertex +1 +1 -1 
vertex -1 +1 +1 
endloop 
endfacet 



endsolid CUBE 
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